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A constant outward radial thrust acceleration can be used to reduce the 
radius of a circular orbit of specified period. Heliocentric circular 
orbits are designed to match the orbital period of Earth or Mars for 
various radial thrust accelerations and are defined as synchronous 
orbits. Minimum-time solar sail orbit transfers to these synchronous 
heliocentric orbits are presented. 


INTRODUCTION 


The concept of using a solar sail as a means of propulsion for a space vehicle was first introduced 
in 1924 by two Russian space pioneers, Konstantin Tsiolkovsky and Fridrickh Tsander. In their theory, the 
transfer of momentum that occurs when a photon of light bounces off the reflective surface of a solar sail 
should provide a mechanism for space propulsion. The effect of this sunlight pressure on the orbital 
motion of interplanetary spacecraft was first taken into account in 1974 with the Mariner 10 missions to 
Mercury and Venus. Solar panels on the Mariner 10 spacecraft when tilted at various angles to the sun 
experienced differences in sunlight pressure, thus inducing rotational motion of the spacecraft. 

In the late 1970’s, the NASA Jet Propulsion Laboratory (JPL) initiated a team effort to study the 
feasibility of using a solar sail vehicle to rendezvous with Halley's Comet during its 1986 approach to the 
inner solar system. Due to the large energy requirements of the mission, both solar sails and solar electric 
propulsion were studied as feasible spacecraft propulsion techniques. Solar electric propulsion was 
ultimately chosen for the mission since solar sails were deemed to not be a sufficiently mature technology 
for a near term space mission. In the end, the U.S. Halley mission was cancelled and research efforts 
directed at producing an operational solar sail vehicle for NASA were halted in 1981. 1 

In an era of tight space mission budgets, solar sail spacecraft fit NASA’s “better, faster, cheaper” 
paradigm and look even more promising than they did in the 1970’s. New materials are accessible today 
for producing extremely thin, strong sheets of sail material. The booms and struts that stretch out and 
maintain a sail’s shape could be made from composite material that is currently available. Also, 
microsatellite technology allows the solar sail payload to have a smaller mass. Together, these new 
technologies add up to feasible solar sail missions. 2 
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This paper considers the mission of transferring between a circular orbit of one astronomical unit 
(AU) and a circular synchronous orbit in minimum time using solar sail propulsion. A synchronous orbit is 
an orbit that matches the period of a celestial body but at a different orbit radius than would exist for an 
non-thrusting spacecraft in a two body system. The ability of a solar sail to provide a constant radial thrust 
allows the existence of these synchronous orbits and will be described in more detail in a later section. 3 
Several investigations into solar sail optimal flight paths have been made in the past. Tsu 4 and London 5 
performed pioneering work using logarithmic spirals to approximate the motion of the sail. Kelley 6 , 
Cavoti 7 , Zhukov and Lebedev 8 and Jayaraman 9 have formulated interplanetary minimum-time transfers for 
solar sail spacecraft using a variety of optimization approaches. Following the methodology of Zhukov and 
Lebedev, Sauer 10 used calculus-of-variations to formulate the necessary conditions for interplanetary 
transfers. Sauer however modeled three-dimensional transfers and included inclined, elliptic orbits for the 
launch and target planets. In this paper, calculus-of-variations is used to derive the first-order necessary 
conditions for a locally optimal minimum-time transfer using solar sail propulsion. The transfers are 
heliocentric beginning in a circular orbit of one AU and terminating in circular orbits that have the same 
period of the Earth or Mars. The transfers that will be presented are planar but the analysis could be 
generalized to three-dimensional transfers if desired. 


SOLAR SAIL EQUATIONS OF MOTION 

To simplify the analysis, two-dimensional planar transfers were modeled using the polar 
coordinate system shown in Figure 1. In the equations of motion, the symbol r denotes the distance from 

the sun to the solar sail center of mass, <f) is the transfer angle between a reference axis and the sun-sail 

line, u is the radial velocity, v is the tangential velocity, and 9 is the sail control angle between the sun- 
sail line and the sail normal measured positive in the counter-clockwise direction from the sun-sail line. 

The force F generated by the sail depends on the total surface area S f the solar intensity W at one 

AU, the speed of light c, r, and the sail control angle 0 . n The acceleration is given below where m is the 
spacecraft mass. 

Fjm = 2WS cos 2 0/cr 2 = £cos 2 0 ! r 2 (I) 

As a further simplification, the constants W t S , and c have been combined into the characteristic 
acceleration constant denoted by k -2WS/c . Equation (1) shows that the constant k is the maximum 
acceleration that the sail can produce at a reference distance unit of 1 AU. 

Sail Normal 
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The variable X is used for the state vector and is defined in Eq. (2). The equations of motion for 
the solar sail in a two-body heliocentric system can then be written as shown in Eq. (3). 



( 2 ) 
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In this paper, heliocentric transfers are modeled and canonical units used, 

therefore = 1 AU^/TU 2 . Equations (2) and (3) are then combined to obtain the state space 
formulation of the equations of motion. 

= * 3 

x 2 = x 4 /x l 

i 3 = -l/jc 2 + k cos 3 0/x 2 + x 4 fx x ^ ^ 

x 4 = -x 3 x 4 /xj +k sin0 cos 2 0/x 2 

OPTIMAL CONTROL PROBLEM 

Optimal control theory 12 will be introduced and then applied to the problem of minimum-time 
transfer to a specified orbit for a solar sail. The nonlinear differential equation describing the motion of the 
solar sail (Eq. (4)) can be written compactly as shown in Eq. (5). 

X = 7[?(0.e" (0, '] x(t 0 ) given, t 0 < t < t f (5) 

In Eq. (5), the state parameters x(/) are determined by the control parameters, 0 (/) where the dot specifies 
differentiation with respect to time, t . The general optimal control problem is then to determine the 
controls 0 ( t ) that will minimize a performance index J that is specific to the problem being solved. 
The performance index J is a scalar representation of the dynamic system's cost. Cost for achieving 
spacecraft orbits using impulsive propulsion systems can typically be defined by the amount of velocity 
change Avor equivalently by the amount of fuel expended. For the case of the solar sail in which 
continuous thrust with no consumption of propellant is possible, the cost being minimized is the time of 
flight. 


A general performance index Jean be written in the form; 

V 

j=j] [4/f'/] + J/(x(/),e (0 ,t)^ 


(6) 
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In the above cost function, rj represents the scalar terminal cost and / represents the integral cost 
accumulated over a particular time interval. For the transfer of a solar sail to a specified orbit in minimum 
time, / = 0 , T] = t/y and therefore the cost function being minimized is the time of flight t j . An 
equivalent formulation is to define 1=1 and rj = 0 . The Hamiltonian can then be written as; 

n = I + X r x=X T i ( 7 ) 

Here the variable X refers to the vector of Lagrange multipliers, also known as the vector of adjoint 
variables. The Hamiltonian for this problem, using Eqs. (4a-d) and four Lagrange multipliers is; 

H = X { (u)+X 2 (v / r ) + X 3 (v 2 / r -1 / r 2 + k cos 3 0 /r 2 )+X 4 (-wv/r + k sin0 cos 2 0 / r 2 ) (8) 

Using the calculus-of-variations optimization technique, a set of necessary conditions can be 
determined. The solution to these necessary conditions provides the 0 (t) that produces a stationary value 
of the performance index J 12 . These first-order necessary conditions are the Euler-Lagrange equations. The 
first of these is the already stated differential equations of motion with the associated boundary condition 
from Eq. (5). The second is the differential equation governing the Lagrange multipliers with the 
associated boundary condition at the final time. These are 

X T =-dH Jdx 

=[dx\ldx+ v T d'F/dl)| ; (9) 

In Eq. (9) X T denotes that the Lagrange multiplier is evaluated at the final time, as is the right hand 
1 V 

side of (9). Also in Eq. (9), 'P is the vector of specified terminal constraints and v~ is an additional set 
of constant multipliers. For transfers to circular synchronous orbits, the three chosen terminal constraints 
v F 1 ,'F 2 , v F 3 are shown in Eq. (10). The values for r specified and v specified will be defined in later sections to 
correspond to desired terminal circular orbits. 


T' = 


^ f specified 

U j - 0 


V / ^ specified 


( 10 ) 


Since there are four states in the X vector, there are four differential equations for X according to 
Eqs. (8) and (9). 

X , = ~dH /dr = —2 A. 3 / /■ 3 +X 3 v 2 //' 2 + 2 kX 3 cos 3 0 /r 3 - X 4 uvfr 2 +2 kX 4 sin0 cos 2 0/r 3 +A 2 v/r 2 

(11a) 

X i(//)=dr| /dry +v ) 3 v F, / dx { +v 2 3 v P 2 Id r f +v 3 S v F 3 /dr 3 = v j =constant (1 lb) 

X 2 =r5H/5* =0 (12a) 


4 


+v,5 v F, /dtj +v 2 0 l F 2 IStj +v i d'¥ i /d$ / = 0 (12b) 

By integration, it is seen that A, 2 ~ constant and by enforcing the boundary condition X 2 = 0 for all 
time. This Lagrange multiplier can be eliminated from the problem formulation if desired leading to a 


reduced order problem. The third and fourth differential equations for X are 

X 3 = -d H/du = X 4 v/r-X l (13a) 

X 3 (/y )= dr\ Iduj -hv , 5 V P 1 Iduj +v 2 5 v T 2 /5wy +v 3 3T' 3 /S*/y =v 2 = constant (13b) 

X 4 = -d H/ d v =-2 \ z v! r + \ 4 u / r (14a) 

X 4 {tf)=&r\/dvj +Vj5 idv j +v 2 d 4^ Idvj +v 3 54 / 3 1 dv j =v 3 = constant (14b) 


The third Euler-Lagrange condition is that the terminal boundary conditions specified in Eq. (10) must be 
satisfied. Since the final transfer time is unknown and free to be determined in a time open problem, the 
fourth Euler-Lagrange equation, the so called transversal ity condition, applies 

[<3r| / dt + //]|^ =0 (15) 

In a time open formulation, dr\ /dt f = 1 , which leads to the transversal ity condition on the Hamiltonian at 
the final time of H(tj ) + 1 = 0 . 


The optimality condition provides a formula for the one control 0 . 

3H/d0=O (16) 

Applying Eq. (16) to the Hamiltonian given in Eq. (8) yields; 

-3X 3 £cos 2 0 sin0 fr 2 +X 4 £(-2cos0 sin 2 0 +cos 3 0)/r 2 =0 (17) 


After dividing Eq. (17) through by £cos 3 0/r 2 a quadratic in tan0 results. 

- 2X 4 tan 2 0 - 3X 3 tan 0 + X 4 = 0 

Equation ( 1 8) can be solved using the quadratic formula to yield the condition for the optimal control 

0 = tan -1 {+[(9X 1+ SX \) V2 -3 X 3 ]/4>. 4 } 


(18) 


(19) 


The control angle in Eq. (19) is uniquely determined by X 3 and X 4 . The sign ambiguity of the square 
root term in the numerator was determined to be strictly positive by examining all possibilities of the signs 
of X 3 and X 4 along with the resultant quadrant where 0 must reside. 13 
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Summarizing, a solar sail minimum-time orbit transfer can be formulated as a two-point 
boundary-value problem defined by the eight differential equations and associated boundary conditions 
Eqs. (4) and (11) through (14) along with the control in Eq. (19). A FORTRAN algorithm was constructed 
to solve this two point boundary value problem. Results obtained from this software were compared with 
published results as described in the next section. 

TEST CASES 

Mars 


In order to validate the computational methodology used to solve the solar sail two point boundary 
value problem, a minimum time solar sail transfer from Earth orbit, r = 1.0 AU, to a Martian orbit, r = 
1.524 AU that is described in Reference 8 was generated. The vector of terminal constraints shown in Eq. 
(10) was driven to zero (within a prescribed tolerance) where /y , Uj , and Vj are the Final values of the 

radius, radial velocity, and tangential velocity and r spec and v spec are the specified values of radius and 
tangential velocity for the desired orbit. For the transfer to a Martian orbit, the values r spec =1 .524 AU and 

v =1 / Vi. 524 were used. A characteristic acceleration of k = 2 mm/s 2 was used to correspond with 

the results given in Reference 8. 

A multi-dimensional shooting method was used to solve the two-point boundary value problem. 
Initial guesses for the transfer time and the Lagrarige multipliers are required with this technique. The 
values of a final time of .5 TU and X ] =1, X 3 =I,andX 4 =1 were used. In the previous section, the 
variable X 2 was shown to be zero for all time. By reinserting successive values of the initial values for the 
three Lagrange multipliers, a converged solution was obtained. This method is known as the “method of 
continuation”. In the converged solution, the constraint vector ¥ achieved the final desired values for r, 
u, and v within a prescribed tolerance of 1 x 10~ 6 . This meant that the final radial velocity u was within 
lxltT 6 of its desired value of zero and that the final tangential velocity v was within lxlO" 6 of its desired 

value of 1 / Vl .524 . Also, the final radius r was within I xl(T 6 of its desired value of 1.524 and thus within 
the Mars sphere of influence (SOI) of 170 Mars radii. The accuracy of the final value of the transversal ity 

constraint was also on the order of lxlO -6 as compared to its ideal value of zero. The final transfer angle 
value was allowed to vary freely. A total minimum transfer time of 324 days or 5.57 TU was obtained fora 
converged solution, which compared favorably with the value of 5.54 TU from Reference 8 for transfer at a 

characteristic acceleration of 2 mm/s 2 . Figure 2 illustrates the control angle state history. Figure 3 
illustrates the control angle at each time interval on the resultant outward spiral trajectory. The state 
histories obtained essentially duplicated the solar sail Mars transfer state histories of Figures 2 and 3 of 
Reference 8. 
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Figure 3. Control Angle Theta at Time Intervals for Solar Sail Earth to Mars Minimum Time Transfer 


Venus 


For a time open transfer to Venus orbit, r =0.7233 AU , Uj = 0 , and v spec = 1 / Vo.7233 were 

used in the terminal constraint vector. A characteristic acceleration of k = 2 mm Is 2 was again used for 
comparison with results from Reference 8. Lagrange multipliers corresponding to a converged solution 
were obtained once again using the method of continuation. In the converged solution, the constraint 

vector ¥ achieved the final desired values for r, u, and v to within a prescribed tolerance of lx 10~ 6 in a 
manner similar to'how the constraints were satisfied for the Martian transfer described above. A total 
minimum transfer time of 163.65 days or 2.815 TU was obtained for a converged solution, which again 
compared very favorably with the value of 164 days from Reference 8 for transfer at a characteristic 

acceleration of 2 mm! s 2 . Figure 4 illustrates the control angle state history. Figure 5 illustrates the 
control angle at each time interval along the resultant inward spiral trajectory. The state histories obtained 
essentially duplicated the solar sail Venus transfer state histories of Figures 4 and 5 of Reference 8. 


The optimization software was able to achieve converged solutions for both the Earth-to-Mars and 
Earth-to Venus solar sail transfers with transfer times and state histories that duplicated the results of 
Reference 8 , the software was considered^ validated. Next, the software was used to investigate inward or 
outward minimum time spirals from EartfTto synchronous orbits of any specified radius. 
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Figure 5. Control Angle Theta at Time Intervals for Solar Sail Earth to Venus Minimum Time Transfer 
SYNCHRONOUS CIRCULAR ORBITS 

The boundary conditions on the radius and velocity magnitudes for synchronous orbits will be 
derived in this section. The solar sail mission statement is to determine the control profile to transfer from 
a circular orbit of 1 AU to a synchronous orbit of a given period in minimum time. The sail control 
authority is governed by the solar sail constant k = IWSjc . Using Eq. (I), k can be shown to be equivalent 
to the maximum acceleration that the sail can provide at 1 AU. To design a synchronous orbit, a desired 
period Paired for a circular orbit is chosen and is written in terms of the radius (r syn ) and velocity (v syn ) of 
the circular synchronous orbit. 

P desired ~ '^^synj^syn ( 20 ) 

The period P desired is next related to the earth orbital period by the factor R. 

P = P desired / P i lutrth (2 1 ) 

Using canonical units with distance referenced to 1 AU, the period of the earth is P carnh = 2 k TU 
and Eq. (20) and (21) are combined to yield a simple relationship between r syn and v syn . 

V V« =r X} ,,/R ( 22 ) 
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Next, the constant radial acceleration that the solar sail provides at the synchronous distance must 
be related with the synchronous velocity. This is done through the equations of motion by enforcing zero 
radial velocity and acceleration on the synchronous orbit. 

v l„ = (m -*)/r 9W = 0 (23) 

Note that by the spacecraft thrusting in only in the outward radial direction once in the 
synchronous orbit, the sail is in essence reducing the local gravitational constant, p. Finally, by combining 
Eqs. (22) and (23), an expression relating the maximum sail acceleration at one AU (k) and the desired 
period is determined. 


r sy „=R 2/3 (]-k) ]n (24) 

The radius of the synchronous orbit supplied by Eq. (24) is used for r speciflcd in Eq. (10). The 
velocity is specified by Eq. (22) and is set equal to v 5p€C }fi C d in Eq. (10). Equations (22) and (24) serve as 
terminal boundary conditions for the transfers along with the sail having zero radial velocity at the end of 
the transfer (i/=0). 


TRANSFERS TO EARTH SYNCHRONOUS ORBITS 

Figure 6 is a plot of Eq. (24) for Earth synchronous orbits, i.e. R=l. Values of k between 0-0.34 
AU/TU 2 (0-2.0162 mm/s 2 ) are presented. This range of sail accelerations was selected due to the fact that 
they are believed achievable using current technology. Figure 6 shows that as the thrust acceleration is 
increased, the synchronous orbit radius moves closer to the sun. For Earth synchronous orbits v Vrt =r J>w . 
The transfer times to the earth synchronous orbits decrease as the thrust acceleration increases even though 
the radial difference between the initial and terminal orbits increases. Points A and B have been labeled on 
Figure 6 corresponding to an acceleration of 0.006 mm/s 2 and 2 mm/s 2 respectively. Transfers to Earth 
synchronous orbits at these two acceleration levels are presented in Figure 7. Converged solutions for 
thrust levels between points A and B were obtained with a tolerance of 10* 6 as in the test cases. 
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k (mm/s 2 ) 

Figure 6 Earth Synchronous Radius and Transfer Times Versus Sail Acceleration 


The transfer to an Earth synchronous orbit corresponding to a 0.006 mm/s 2 is barely 
distinguishable from the circular reference orbit at 1 AU. The transfer requires 182 days, approximately 
one-half of the period of the circular reference orbit. The transfer at the higher acceleration (2 mm/s 2 ) 
requires only 152 days to complete the transfer. 


The control profiles for these two transfers are displayed in Figure 8. Each transfer shows similar 
characteristics. A large slew maneuver is required to match the desired boundary conditions on the 
velocity. The maneuver for the lower acceleration occurs at approximately the middle of the transfer while 
the maneuver in the larger acceleration transfer is slightly delayed. Note that the control angle is always 
negative corresponding to the acceleration being applied in the direction opposite to the velocity. This 
control decreases the orbital energy of the spacecraft causing an inward spiral to occur. 


12 










TRANSFERS TO MARS SYNCHRONOUS ORBITS 


Figure 9 is a plot of Eq. (24) for Mars synchronous orbits, i.e. R=(ly 321 .73d)/(ly)=l. 88083. 14 
Figure 9 shows that as the thrust acceleration is increased, the synchronous orbit radius moves closer to the 
Earth. Also, recall that for synchronous orbits v, >7 ,=r> >7 /R. The transfer times to the Mars synchronous 
orbits rapidly decrease as the thrust acceleration increases (See Figure 10) since the radial difference 
between the initial and terminal orbits decreases along with the increased control authority. Points C and D 
have been labeled on Figure 9 corresponding to an acceleration of 0.38 mm/s 2 and 2 mm/s 2 respectively. 
Converged solutions for thrust levels between points C and D were obtained with a tolerance of I O' 6 as in 
the test cases. Minimum-time transfers for characteristic accelerations less than 0.38 mm/s 2 were not 
obtained due to the strict convergence requirements and the fact that the sail has less control authority and 
larger energy changes than the higher acceleration cases. f 
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Figure 9 Mars Synchronous Radius and Velocity Versus Sail Acceleration 
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Figure 10 Mars Synchronous Transfer Times Versus Sail Acceleration 


Figure 1 1 displays the trajectories corresponding to points C and D. The transfer to an Mars 
synchronous orbit corresponding to a 0.38 mm/s 2 travels a transfer angle of 320 degrees and requires 429 
days complete. The transfer at the higher acceleration (2 mm/s 2 ) requires only 123 days to complete the 
transfer. 


The control profiles for these two transfers are displayed in Figure 12. Due to the small 
acceleration of point C, the angle between the sun-sail line and the sail normal is always positive (thrusting 
in the positive velocity direction ) in order to increase the energy. The control angle history for point D is 
quite different. Due to the larger available control, the boundary conditions are reached in minimum time 
by thrusting away from the velocity direction. 

CONCLUSIONS 

Minimum time transfers for solar sail spacecraft to Earth and Mars synchronous heliocentric orbits 
have been determined for sail accelerations that are achievable using current technology. Necessary 
conditions for local optimality were derived using calculus-of-variations and boundary conditions for 
synchronous circular orbits were presented. A multi-dimensional shooting method was implemented to 
solve the boundary value problem. The optimization algorithm was verified on two published 
interplanetary transfers and then applied to the minimum-time synchronous transfers. Solar sail spacecraft 
in heliocentric Earth synchronous orbits could be used to warn spacecraft of solar flare events that could 
have significant effects for power generation and communication on Earth. Mars synchronous orbits would 
provide locations for communication satellites with constant Earth-Mars viewing angles. Solar sail 
spacecraft have a great deal of potential for advancing the state-of-the-art in spacecraft propulsion 
technologies and the further exploration of the solar system. 
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Figure 1 1 Transfers to Two Mars Synchronous Orbits 
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Figure 12 Control Angle to Two Mars Synchronous Orbits 
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